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Abstract. We investigate the dynamics of and the relations between small-scale penumbral and photospheric 
features near the outer penumbral boundary: penumbral grains (PCs), dark penumbral fibrils, granules, and 
photospheric G— band bright points. The analysis is based on a 2 h time sequence of a sunspot close to disc center, 
taken simultaneously in the G— band and in the blue continuum at 450.7 nm. Observations were performed at 
the Swedish Vacuum Solar Telescope (La Palma) in July 1999. A total of 2564 images (46" x 75") were corrected 
for telescope aberrations and turbulence perturbations by applying the inversion method of phase diversity. Our 
findings can by summarized as follows: (a) One third of the outward- moving PCs pass through the outer penumbral 
boundary and then either continue moving as small bright features or expand and develop into granules, (b) Former 
PCs and G— band bright points next to the spot reveal a different nature. The latter have not been identified 
as a continuation of PCs escaping from the penumbra. The G— band bright points are mostly born close to dark 
penumbral fibrils where the magnetic field is strong, whereas PCs stem from the less-magnetized penumbral 
Q ' component and evolve presumably to non-magnetic granules or small bright features. 
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1. Introduction granules. However, the granular convection is perturbed 

by the presence of many magnetic elements, visible in 

$ The magnetic field generated at the bottom of the con- G _ band as tiny bright points located in Agranular 

vection zone emerges m active regions, from where it is gpaces _ They correspond to the facular and magne tic ele- 

diffused to the surface of the Sun, by convective motions ments which are yisible in lower resolution mtergrams and 

and large scale plasma circulation. However, the mecha- magne tograms. The granules (Muller & Mena 1987; Simon 

msm by which magnetic flux is taken away from sunspots et aL 19g8) and the magnetic e i eme nts - either observed 

is not well known yet. The aim of this work is to investi- in nltergrams or in magnetograms - (Sheeley 1969, 1972; 

gate how the features inside and outside the penumbra of shedey & Bhatnagar 1971; V rabec 1971, 1974; Harvey 

a decaying sunspot are related. Decaying sunspots are well & Haryey lg73; Muller & M&a ig87; Brickhouse & La 

suited for that purpose, because it is believed that the de- Bonte 1988; Lee 1992; R yut0 va et al. 1998 ) move radially 

cay of a sunspot proceeds predominantly by erosion of its outwardS; through an annular cell. This annular cell has 

perimeter leading to a quadratic decay law (Sheeley 1972; been called a « moat » and the magnet ic elements, moving 

Meyer et al. 1974; Petrovay & Moreno-Insertis 1997). acC ross, "moving magnetic features" or "MMFs" ( Sheeley 

Sunspots are very complex magnetic structures, em- 19 7 2; Harvey & Harvey 1973). Dopplergrams have shown 

bedded in a convective plasma, which show a small-scale tbat tms outflow, with speed of 0.5-1.0 km s~\ is very 

pattern, down to the resolution limit of the largest contem- s i mnar to that observed in supergranules. This has sug- 

porary solar telescopes (-O'.'l = 75 km). The surround- ges t e d that the moat is presumably a supergranule whose 

ing convective plasma is mainly observed in the form of center is occupied by a sunspot, and which extends ap- 
proximately 10000 km beyond the spot's edge (Sheeley 
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1972) . The penumbra-granulation boundary appears very 
sharp, since the photospheric features (granules and bright 
points) are directly adjacent to darker penumbral elements 
(dark diffuse patches, dark and bright fibrils). 

The penumbra has a very complex morphological, 
magnetic and dynamic structure, which is not yet properly 
understood physically. Visually, with the spatial resolu- 
tion of the present work (~0'/23), it appears to be formed 
by narrow, nearly radial bright filaments, separated by 
dark fibrils. Local brightenings, called penumbral grains 
(PGs) are observed in the bright filaments. In images of ex- 
tremely high resolution (better than 0'.'12) obtained with 
the 1-m Swedish Solar Telescope at La Palma, Scharmer et 
al. (2002) detected dark cores inside the bright filaments. 

The magnetic field in the penumbra is as strong as 
1500 G at the inner boundary, and decreases to 700 G 
at the outer boundary, becoming more and more hori- 
zontal. In addition, it is very inhomogeneous in the az- 
imuthal direction. From the synthesis of many papers, 
sometimes presenting conflicting results, it appears that, 
in the outer penumbra, higher field strengths tend to be 
correlated with more horizontal fields (by 10 to 20 degrees) 
and dark intensity structures (see e.g. Beckers & Schroter 
1969; Westendorp Plaza et al. 2001; Bellot Rubio 2003 and 
references therein). In addition to that, several authors re- 
port that the magnetic radius of sunspots extends beyond 
the visible boundary of the penumbra (Skumanich 1992; 
Title et al. 1993; Martinez Pillet 1997; Westendorp Plaza 
et al. 2001; Bellot Rubio 2003). 

The penumbra of sunspots is not only a magnetic, but 
also a very dynamic plasma: The bright PGs and dark fea- 
tures both move, either toward or away from the umbra; 
an irregular and not steady outflow of gas (the Evershcd 
effect) is detected. In the inner penumbra, the PGs always 
move towards the umbra, and reach a maximum velocity 
of 0.5 km s" 1 at the umbra-penumbra boundary (Muller 

1973) . In the outer penumbra, they have been found to 
move either towards the surrounding photosphere (Wang 
and Zirin 1992; Denker 1998; Sobotka et al. 1999; Sobotka 
& Siittcrlin 2001), or towards the umbra (Muller 1973; 
Tonjes & Wohl 1982; Zirin & Wang 1989). The divid- 
ing line, separating the two opposite motions, is located 
near the outer boundary at about 1/3 of the width of 
the penumbra. Schlichenmaier (2002) elaborated a mov- 
ing flux tube model developing waves and kinks (a pho- 
tospheric "serpent") that reproduces both directions of 
motions of the PGs. 

The present paper is the first in a series based on 
two excellent 2-hour time series of sunspot images recon- 
structed with the phase-diversity technique. The long du- 
ration of these series, together with the high and stable 
quality throughout the whole period, which is substan- 
tially improved after the reconstruction (see Fig. 1), makes 
this material one of the best data sets ever produced to 
study the morphology and the dynamical behaviour of a 
sunspots fine structure and its surroundings. In addition, 
the strict simultaneity of the images in the G— band and 



the blue continuum significantly increases the reliability of 
the identification and tracking of G— band bright points. 

In this work we are interested in the relation between 
the penumbral and the photospheric features as they move 
across the boundary of a decaying sunspot. Such a decay- 
ing sunspot is well suited for our purpose, because it has 
a simple, nearly symmetric configuration, surrounded by 
a well developed and nearly circular moat. Moreover, the 
features in the outer penumbra and in the moat move in 
the same direction, outward, which is a favourable situa- 
tion to find out if there is a continuity in their properties. 
In the second paper of the series, we will investigate in 
detail the motions of granules, families of granules formed 
by recurrently splitting granules, and small-scale magnetic 
bright points, in order to see how they behave and inter- 
act when they are dragged by a supergranular convective 
flow. 

2. Observations and Data Processing 

2.1. Observations 

The leading sunspot of a decaying bipolar group NOAA 
8620 was observed on 7 July 1999 at the Swedish Vacuum 
Solar Telescope, La Palma (Scharmer et al. 1985), near the 
center of the solar disc at /i = 0.93. The sunspot showed 
a regular shape and a slow and smooth area decay from 
July 3 to July 13 as recorded in the daily area values of 
Solar Geophysical Data. 

Two strictly simultaneous time sequences of high- 
resolution images were taken in the wavelength bands at 
450.7 ± 0.45 nm (blue) and 430.8 ± 0.55 nm (G-band). 
The observing sequence spans for more than two hours 
(1282 images per series), from 7:30 UT to 9:43 UT. The 
images were acquired using two Kodak Megaplus 1.6 
CCD cameras, with a dynamical range of 10 bit and 
1536 x 1032 pixels. The pixel size and exposure time were 
0'.'083 and ^40 ms, respectively. The cameras were placed 
in orthogonal light paths produced by a cubic beam split- 
ter. A common mechanical shutter in front of the beam 
splitter ensured the simultaneous exposure in both cam- 
eras. 

To apply the phase-diversity technique for image re- 
construction, pairs of focused and intentionally defocuscd 
images must be taken simultaneously To that aim each 
CCD camera is equipped with a special beam-splitter that 
images the focused-defocused image pair onto the same 
chip (see Fig. 2 in Lofdahl et al. 1998). The defocus is 
caused by a slightly longer optical path length in one of the 
split beams. In the present case, the equivalent defocusing 
displacement in vacuum was 8.85 mm that corresponds to 
a phase shift at the edge of the aperture of 1.15 waves and 
l.lOwaves at 430.8 nm and 450. 7nm, respectively. After 
selection of the common area in both images of the pair, 
the removal of spurious effects induced by the beam split- 
ter, and by the apodization during restoration process at 
the edges of the image, the effective field-of-view results 
in 46" x 75" (560 x 910 pixels). The cameras were work- 
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ing in a real-time frame-selection mode storing the best 
4 images within a time interval of 18 s. Considering the 
additional time to transfer the data to the hard disc, the 
resulting mean time interval between two consecutive sets 
of 4 images was 25 s. 

In an independent third channel, a complementary 
data set was acquired in the Ca II K line (A 393.3 ± 0.15 
nm ) for monitoring the magnetic activity in the observed 
field. Standard flat-field and dark-current corrections were 
applied to all images before the image restoration. 

2.2. Data Processing 

As mentioned above, the phase-diversity (PD) technique 
has been applied to correct the images of the blue and 
G— band sequence up to close the diffraction limit. This 
technique allows one to estimate both the unknown ob- 
served object and the wavefront aberration. A few pa- 
pers in the literature are complementary in the descrip- 
tion of the PD technique for image reconstruction (e.g. 
Gonsalves & Childlaw 1979; Gonsalves 1982; Paxman, 
Schulz & Fienup 1992; Lofdahl & Scharmer 1994: and 
Paxman et al. 1996). However, we summarize the most 
relevant items of the method as an introduction to this 
powerful tool for high resolution solar imaging and to re- 
port on our particular setting of parameters and assump- 
tions to run the computing code. 

This technique requires the detection of at least two 
images of the object. One of these is the conventional 
focal-plane image that is degraded by unknown aberra- 
tions induced by the atmosphere and the telescope. The 
other one is a simultaneous image affected by the same un- 
known aberrations plus an extra known aberration (the 
easiest controlled aberration to produce is a defocus). 
Both image realizations, ii(q) and i 2 (q), can be math- 
ematically represented by the system of equations 



n(<z) 



«o(<?) * si(q) 
io(q) * s 2 (q) 



n 2 \ 



(1) 



where * stands for convolution and q is the vectorial 
notation for the coordinates of the image points; io is the 
"true" object and i\, i 2 are the two observed images; s\ 
and s 2 represent the corresponding point spread functions 
(PSF) of the optical transmission system - including the 
terrestrial atmosphere and the telescope - and n\, n 2 are 
additive terms of noise. 

Gonsalves & Childlaw (1979) propose the solution of 
this system of equations as a least-squares fit which in 
the Fourier domain can be written as the following error 
metric: 



u l 

+ 1 \i 2 ( U )-i (u)s 2 (u)\ 2 y 



(2) 



where capital letters represent the Fourier transform of 
the corresponding functions denoted by lower-case charac- 
ters and u is the frequency vector in the Fourier domain. 7 



is included as a factor to equalize the noise in both images 
(a ratio of the respective noise variances). Typical values 
for noise level in our data at both wavelengths are in the 
range 0.6-0.8% of the average intensity. S represents the 
optical transfer function (OTF). It is directly related to 
the phase aberration in the pupil (see e.g. Bonet 1999), 
which in turn can be approximated by the first J terms 
of a Zernike expansion. In the present case Zernike modes 
2-21 have been employed in the expansion, i.e. 20 coeffi- 
cients, a = {ctj,j — 2, ...,21} have to be determined to 
estimate the wave front error and subsequently S(u,a). 
Using a greater number of modes, e.g. J = 45, did not 
significantly improve the results while it substantially in- 
creased the computing time. 

Part of the minimization of J5J) can be performed an- 
alytically. The solution of equation: dL/dlo — 0, gives an 
estimate of the object, Iq, that minimizes for fixed a. 



Io(u) 



h(u) S{(u, a) + 7/2(11) Sl(u, a) 
|5 1 (t t ,a)| 2 + 7 |5 2 («,a)| 2 ■ 



(3) 



The subsequent substitution in (J2J leads to the modified 
error metric Lm{ol) 



Lm{o) = 



u 



\h(u) S 2 (u, a) - I 2 (u) Si(u, a) 
|Si(«,a)| 2 + 7 |S 2 (ix,a)| 2 



(4) 



which is now independent of the object Fourier transform 
Iq. Thus, the parameter space over which the optimization 
is performed has dimension 20. Once these parameters are 
determined, Si and S 2 can be constructed and the object 
estimate, Iq, can be derived from J3J, thereby completing 
the image reconstruction process. 

Poor seeing may produce a low signal-to-noise ratio 
and very small and uncertain values of S(u) at certain 
spatial frequencies. At these frequencies, Eq. J2J) produces 
an excessive amplification of the signal, giving rise to ar- 
tifacts that show up as a regular pattern in the restored 
scene. To circumvent this problem, we combine in a sort 
of speckle summation (see e.g. Paxman et al. 1996) the 
results of the inversion of 3 or 4 image pairs taken close in 
time, so that evolutionary aspects in the solar structures 
can be neglected. The lower right panel of Fig. 1 shows 
the central part (41'.'5 x 33'.'2) of the final restored image 
resulting from the combination of three single-image in- 
versions of the blue sequence. These images were recorded 
within a period of 18 s. 

The restored images were de-rotated to compensate 
for rotation of the field of view in the focal plane in- 
duced by the alt-azimuthal configuration of the telescope. 
Furthermore, the images were aligned, destretched (com- 
puter code by Molowny-Horas & Yi, 1994), and filtered 
for p-modes in the k — u space (cut-off phase velocity 
5kms _1 ). Two movies were produced, for the blue con- 
tinuum and the G— band, of 288 frames each, spanning 
over 2 hours (At = 25 s) and covering a field of 46" x 75", 
with the sunspot located in the center. Direct inspection 
of the images reveals structures up to about the Rayleigh 
resolution limit at both working wavelengths (^0"23). 
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Fig. 1. Sunspot images taken from the sequence at 450. 7nm. Lower right panel: Central part of the restored image 
resulting from the combination of three raw images observed close in time and shown in the other three panels of 
the figure. Each raw image has been individually inverted with phase diversity in order to evaluate its corresponding 
instantaneous OTF. 



A feature-tracking technique (Sobotka et al. 1997) was 
used to measure positions, intensities and sizes of PGs 
and of G— band bright points. PGs were tracked in the se- 
ries of blue images. To isolate them from other structures, 
photospheric granules and umbral dots were masked out 
and a segmentation algorithm was applied to each blue 
image, based on the rule that pixels with downward con- 
cavity, representing bright features, are set to 1 and the 
rest to (cf. Sobotka & Siitterlin 2001). The obtained bi- 
nary masks were then multiplied by corresponding blue 
images, so that original intensities of PGs were preserved. 

To isolate G— band bright points, the destretch algo- 
rithm was applied to each blue image to align it exactly 
with the structures in the corresponding G— band frame. 
The granular contrast of aligned blue images was reduced 
by a factor of 0.833 to match the contrast observed in the 
G— band frames. Then the blue frames were subtracted 
from the G— band frames, thus producing intensity differ- 
ence images with the G— band bright points substantially 
enhanced. The segmentation was done by thresholding 
these difference images, obtaining binary masks, and by 
multiplying the masks by the original G— band frames. 

After the segmentation, the features under study (PGs 
and G— band bright points) are formed from the non-zero 
intensity pixels. In the next step, pixels forming a feature 



are labelled by an identification number. Then, the spatial 
coincidences of features in each pair of subsequent images 
are investigated. Two features are identified as predeces- 
sor/successor if they coincide in the coordinates of at least 
one pixel in both frames. The maximum intensity, its co- 
ordinates and the total number of pixels in the feature are 
recorded for each frame. The lifetimes (number of frames) 
and the sizes (number of pixels) of the features can be ob- 
tained directly from the records. Formation, death, split- 
ting and merging of features are taken into account. In 
the case of splitting, the brightest feature is adopted as 
the successor, while if merging occurs, the merged feature 
is defined to be the continuation of the feature with the 
longest record. 

To improve the reliability of the tracking results and 
to avoid small-scale noise in the segmented images, for 
the further analysis we select only features with lifetimes 
and time-averaged sizes above certain minimum values 
(see Sect. 3), where the minimum size is usually derived 
from the resolution limit of the observations. Variable im- 
age quality may cause spurious merging and splitting that 
can make independent features appear related. This would 
lead to large displacements of the maximum intensity po- 
sitions and to unrealistic velocities of proper motions. This 
problem can be partially fixed by selection of features with 
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time-averaged velocities smaller than the cut-off velocity 
used in the k — ui filter. Time-averaged velocities are calcu- 
lated using linear least-squares fits to the positions. More 
additional criteria can be applied to select features with 
specific characteristics. 

Finally, a "visual consistency check" is applied to the 
selected features: Trajectories are checked visually to elim- 
inate discontinuities in position or strong bends/breaks in 
the trajectory, caused by spurious coincidences of differ- 
ent features, and to discard such features that are not the 
subject of the study. 

3. Results 

3.1. Outward-moving penumbral grains 

We have already mentioned in Sect. 1 that in the outer 
1/3 of the penumbra most of the PGs move towards the 
outer penumbral border (P/G boundary). 

A natural question is what happens when they reach 
the boundary between the penumbra and the surrounding 
photospheric granulation (P/G boundary). A visual 
inspection of a movie composed of the blue images shows 
that the outward-moving PGs evolve in three possible 
ways: 

1. They disappear in the penumbra before reaching the 
P/G boundary. 

2. They cross the P/G boundary and continue to move in 
the granulation as small bright features. 

3. After crossing the P/G boundary, they increase their 
size and develop into granules that move away from the 
sunspot. 

This different behaviour of the PGs crossing the P/G 
boundary is quite remarkable and justifies the separate 
study of the features that are observed to expand and 
those that do not expand. 

To study the evolution of PGs quantitatively, the men- 
tioned above feature-tracking technique was applied to a 
field of 12'.'5 x 20'.'8, including a part of the penumbra 
and adjacent granulation. This region is located on the 
left-hand side of the sunspot in Fig. 1 (see also Fig. 2 
showing the contour of the time-averaged outer penum- 
bral border in this region). The criteria to select bright 
features of interest from the raw results of the tracking 
(see Sect. 2) were as follows: (a) Minimum lifetime of 250 s, 
making it possible to pass a distance of 1000 km with the 
(b) maximum allowed average speed of 4kms _1 ; (c) min- 
imum average size of 0'.'28, corresponding to an area of 
9 pixels. To eliminate from the tracking results most of 
the granules and PGs that probably do not cross the P/G 
boundary, an additional criterion (d) of average brightness 
ranging from / = 0.9 to 1.15 (in units of the mean photo- 
spheric intensity) was applied. This criterion was derived 
from the histogram of time-averaged intensities (resulting 
from the tracking), which clearly showed two populations, 
PGs and granules, with peaks at I = 0.92 and 1.16, re- 
spectively, and with a dip at / = 1.02. We can expect that 



the time-averaged intensities of PGs that cross the P/G 
border and penetrate into the granulation correspond to 
the overlap of both populations. This way, 712 bright fea- 
tures were selected. Finally, the visual consistency check 
was applied to these features and all inward moving PGs, 
together with features originating outside the penumbra 
were discarded. 

After the visual consistency check we obtained a sam- 
ple of 126 outward moving PGs, including their tra- 
jectories. The noise in the measured positions was re- 
moved by smoothing the trajectories using cubic splines. 
Instantaneous velocities were calculated as derivatives of 
the smoothed trajectories. Each PG was inspected indi- 
vidually to check its positions in time with respect to 
the P/G boundary. About 2/3 of PGs disappear before 
reaching the P/G boundary, 1/6 cross the boundary and 
continue moving as small bright features with a diameter 
less than 0'.'5, and about 1/6 cross the boundary, expand 
in size, and develop into granules -here we use the term 
granules because the visual inspection of the movie reveals 
that these structures, when present in the photosphere, do 
behave as normal granules. All the features that cross the 
P/G boundary preserve their outward motion. Their tra- 
jectories, smoothed by cubic splines, are plotted in Fig. 2 
together with the time-averaged contour of the P/G bor- 
der. This contour is shown only for illustration, because 
the ragged border of the penumbra evolves in time and 
only individual tracking of PGs gives the information if 
they crossed the P/G boundary or not. The origins of tra- 
jectories are marked by asterisks for small features and 
by squares for granules. Dots denote the following posi- 
tions; in the case of slow motions the dots merge, forming 
solid lines. In general, the birth places of crossing and 
non-crossing PGs are mixed. However, on average, cross- 
ing PGs in our sample originate at the distance of about 
2" from the P/G boundary while the non-crossing ones 
by Of! 9 deeper in the penumbra (the width of the penum- 
bra is approximately 8"-10"). The present sample of PGs 
does not include all PGs in the field of view. 

The total lifetimes (including the penumbral and pho- 
tospheric stages) of the PGs that convert into small bright 
features are in the range of 4.6-41.0 min with the mean 
of 13.9 min. The total lifetimes of the PGs that develop 
into granules are longer: 5.9-120.0 min, with an average of 
36.0 min. For comparison, the non-crossing PGs live 4.2- 
102.5 min, with a mean lifetime of 14.3 min. After crossing 
the penumbral border, small features have lifetimes in the 
range 2.5-22.1 min with the mean of 7.7 min. For gran- 
ules which arise from PGs, the range is 1.7-36.3 min and 
the mean is 15.0 min. Alissandrakis et al. (1987) give a 
mean lifetime of photospheric granules of 16 min. Hence, 
after crossing the border, granules of penumbral origin 
have a mean lifetime comparable with granules formed 
in the quiet photosphere. Small features live shorter than 
"average" granules but longer than small (< 0'.'5) ones, 
with a lifetime that was estimated at 1 3 min (Kawaguchi 
1980; Hirzbcrgcr et al. 1999). 
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Table 1. Summary of PGs near the P/G border 
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Fig. 2. Trajectories of PGs that crossed the P/G bound- 
ary, approximated in the figure by a time-averaged con- 
tour. Origins of trajectories are marked by asterisks for 
small features and by squares for granules. The coordi- 
nate unit is 1 pixel, i.e., 0'.'083. 



The lifetime-averaged speeds of small bright features 
range from 0.4 to 2.7kms~ 1 with a mean of 1.4kms _1 . 
The newly formed granules are slower: 0.2-1.9 kms -1 , 
with a mean speed of 1.1 kms -1 . These velocities are 
slightly higher compared to the average speed of all 
outward moving PGs near the P/G border (0.9 kms -1 , 
Sobotka & Suttcrlin 2001) and to the horizontal speed 
of granules close to sunspots (0.5-0.8 kms -1 , Muller & 
Mena 1987, Simon et al. 1988). 

Figure 3 demonstrates the transformation of an out- 
ward moving PG into a small bright photospheric feature 
(top) and into a granule (bottom). Dark dots in the figure 
denote intensity maxima of the features. 

In Table 1 we summarize the numbers of PGs that dis- 
appear in the penumbra before crossing the P/G bound- 
ary, PGs that cross the P/G border and convert into small 
bright photospheric features, and those that develop into 
granules. We also show the numbers of PGs that accelerate 
or decelerate their motion before and during the crossing 
of the boundary. Note that most of the PGs which con- 
vert into small features decelerate before they escape from 
the penumbra, while those that expand and develop into 
granules accelerate. 



PGs 


Total number 


Accel. 


Decel. 


Non-crossing 


85 (67%) 






— > small features 


21 (17%) 


6 


13 


— > granules 


20 (16%) 


11 


6 



3.2. Relation of G— band bright points to the 
penumbra 

Numerous G— band bright points were observed in the 
sunspot moat, which, in our case, had approximately an 
annular shape and an average width of 8". This has been 
determined from the horizontal velocity field of granules 
around the spot. The details will be described in the sec- 
ond paper of this series. The feature tracking algorithm 
was applied to the series of segmented G— band frames, 
detecting 776 G— band bright points with lifetimes longer 
than 2.5 min and with time-averaged sizes larger than 
0'.'28.Most of the G— band bright points (95 %) move away 
from the sunspot. A detailed analysis of their trajectories 
and horizontal velocities will be published in the second 
paper of this series. 

To find a possible relation of the G— band bright points 
with the penumbra and penumbral features, we looked for 
such G— band bright points that were born in the neigh- 
bourhood of the penumbral border. The neighbourhood 
was defined by the maximum distance of 0'.'25 (3 pixels) 
between the birth position of the G— band bright point and 
the nearest penumbral feature, i.e, PG or a dark penum- 
bral area or a dark fibril. The distance of 3 pixels is com- 
parable with the Rayleigh's resolution limit (0'.'23) of a 
diffraction-limited 47.5 cm telescope at 430 nm. The birth 
positions of 776 G— band bright points, obtained from the 
feature tracking, were compared visually (see examples 
in Fig. 4) with the locations of penumbral features in 
the aligned blue images (Sect. 2). The blue wavelength 
band was selected because of higher contrast. In total, we 
found 191 G— band bright points born in the neighbour- 
hood with the penumbra and then drifting away from the 
spot. Most of them, 132 (69%), originate in the vicinity of 
dark penumbral areas or fibrils, where the magnetic field is 
expected to be stronger and more horizontal than in bright 
filaments. Very often, the G— band bright points are born 
on the tips (ends) of dark penumbral fibrils. Thus, in the 
vicinity of a decaying sunspot, most of the G-Band bright 
points, which are believed to be associated with thin mag- 
netic flux tubes, appear to be born close to a magnetized 
plasma. Only 59 (31 %) G— band bright points appear near 
bright penumbral features (PGs or bright filaments) or 
close to diffuse regions in the penumbra, which are dif- 
ficult to classify as "dark" or "bright". Fig. 4 represents 
blue continuum images showing the birth positions of six 
G— band bright points. Since they are not G— band images 
the bright points are not evident in the figure and, con- 
sequently, white circles -with the radius of 0'.'25- are em- 
ployed to indicate their initial positions or equivalently to 
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Fig. 3. Series of frames showing two examples of PGs crossing the P/G boundary and converting into a small bright 
feature (top) and into a granule (bottom). The tracked features are marked by small black dots at the positions of 
intensity maxima. Numbers in the upper left corners show the elapsed time in seconds. The field of view is 12'.'5 x 
7"5. 



outline their neighbourhood. Black lines represent trajec- 
tories of the bright points' motions. Four of these G— band 
bright points were born on the tips of dark penumbral fib- 
rils and two in the vicinity of bright PGs. 

The small bright features escaping from the penum- 
bra and the G— band bright points born close to the P/G 
boundary are, at first sight, similar in morphology and dy- 
namics. To check for possible coincidences between them, 
we compared, frame by frame, the positions of PGs after 
crossing the P/G border with the positions of all G— band 
bright points. We did not find any case of coincidence in 
time and position, so that we can conclude that escaped 
PGs and G— band bright points are different types of fea- 
tures and that G— band bright points that are located near 
the P/G boundary do not originate from bright features 
which crossed the P/G boundary. 

In quiet granulation, G— band bright points usually 
appear in dark spaces which are compressed by ex- 
pansion or convergent motions of neighbouring granules 
(Muller 1983; Mullcr ct al. 1989; Muller & Roudier 1992). 
Inspecting the G— band movie, we observe that in 95 % 
of all cases, the formation of G— band bright points in the 
neighbourhood of the penumbral border also occurs in ar- 
eas compressed by convergent horizontal motions caused 
by outward moving penumbral features and expanding 
granules adjacent to the border. In quiet granulation, the 
formation of G— band bright points can be attributed to 
the concentration of diffuse magnetic flux, already present 
in the intergranular space, by converging or expanding 
adjacent granules. A similar mechanism of concentration 
of magnetic flux could also take place in the magnetic 
penumbral plasma at the P/G border. 



Fig. 4. Blue continuum images showing the birth positions 
of four G— band bright points (top) born near the tips of 
dark penumbral fibrils and two (bottom) in the vicinity 
of bright PGs. G— band bright points arc not evident in 
the blue continuum so that they are represented by white 
circles with radii 0'.'25 marking their initial positions and 
by black lines representing their trajectories. The field of 
view in the images is 12'.' 5 x 8'.' 3. 

4. Discussion and conclusions 

Two time series of 288 frames each, spanning over two 
hours and taken simultaneously in the G— band and the 
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blue continuum were restored for instrumental and atmo- 
spheric degradation using the phase-diversity technique. 
The excellent quality of the restored images and the stabil- 
ity with time made it possible to study relations between 
small-scale penumbral and photospheric features near the 
outer penumbral boundary: PGs, dark penumbral fibrils, 
granules, and G— band bright points. 

It was shown that about 1/3 of outward moving PGs 
located near the outer penumbral boundary escape from 
the penumbra and penetrate into the surrounding granu- 
lation where they continue their outward motion, either as 
small bright features, or growing as expanding granules. 
Their velocities are slightly higher than average speeds in 
the outer penumbra and in the sunspot moat. They do 
not evolve as G— band bright points in the sense that the 
PGs escaping from the penumbra never coincide spatially 
with the G— band bright points that are located near the 
P/G boundary. 

While some PGs, belonging to the less-magnetized 
component in the penumbra, convert into presumably 
non-magnetic or weakly magnetized granules, many 
G— band bright points are born close to dark penumbral 
fibrils (or directly at their tips), where the magnetic field 
is strong and nearly horizontal. This confirms the finding 
by Title et al. (1995) and Shine et al. (1996), reported 
by Ryutova et al. (1998), that "moving magnetic fea- 
tures" frequently appear along the continuation of dark 
filaments. In addition, our observations show that local 
convergent motions, caused by the expansion of granules 
and outward moving penumbral features, may play a role 
in their formation. The compression of a dark space by the 
surrounding granules, either in the photosphere away from 
sunspots, or at the sunspot border, should be accompa- 
nied by a strong downflow, which may help to concentrate 
the flux into a magnetic flux tube, detected as a moving 
magnetic feature or a G— band bright point. It would be 
of much interest to obtain high spatial resolution spectra 
simultaneously with G— band images, to check whether 
the formation of G— band bright points is associated with 
strong downflows in dark filaments at the P/G border of 
sunspots. 

We have not found in the literature any sunspot 
penumbra model able to explain our observations. 
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